Quantum Fluctuations, Temperature and Detuning Effects in Solid-Light Systems 
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The superfluid to Mott insulator transition in cavity polariton arrays is analyzed using the variational cluster 
approach, taking into account quantum fluctuations exactly on finite length scales. Phase diagrams in one 
and two dimensions exhibit important non-mean-field features. Single-particle excitation spectra in the Mott 
phase are dominated by particle and hole bands separated by a Mott gap. In contrast to Bose-Hubbard models, 
detuning allows for changing the nature of the bosonic particles from quasilocalized excitons to polaritons to 
weakly interacting photons. The Mott state with density one exists up to temperatures T/g > 0.03, implying 
experimentally accessible temperatures for realistic cavity couplings g. 
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The prospect of realizing a tunable, strongly correlated sys- 
tem of photons is exciting, both as a testbed for quantum 
many-body dynamics and for the potential of quantum sim- 
ulators and other advanced quantum devices. Three propos- 
als based on cavity-QED arrays have recently shown how this 
might be accomplished til |2[ [3D, followed by further work 
Engineered strong photon-photon interactions 
and hopping between cavities allow photons (as a component 
of cavity polaritons) to behave much like electrons or atoms 
in a many-body context. It is clear that a particular signa- 
ture of quantum many-body physics, the superfluid (SF) to 
Mott insulator (MI) transition, should be reproducible in such 
systems and be similar to the widely studied Bose-Hubbard 
model (BHM). Yet, the BH analogy is not complete. The 
mixed matter-light nature of the system brings new physics 
yet to be fully explored. 



"Solid-light" systems — so-named for the intriguing MI 
state of photons they exhibit — are reminiscent of cold atom 
optical lattices (CAOL) |9|], but have some advantages con- 
cerning direct addressing of individual sites and device inte- 
gration, and the potential for asymmetry construction by in- 
dividual tuning, local variation, and far from equilibrium de- 
vices [4]. Photons as part of the system serve as excellent 
experimental probes, and have excellent "flying" potential so 
that they can be transported over long distances. Tempo- 
ral and spatial correlation functions are accessible, and non- 
equilibrium quantum dynamics may be studied using coherent 
laser pumping to create initial states. The possible implemen- 
tations are man y In particular, microcavities linked by 
optical fibers lElllOll. small arrays of stripline superconducting 
Cooper-pair boxes or "transmon" cavities OilLUJ], condensate 
arrays 111211 and color center/quantum dot periodic band-gap 
(PBG) materials [01 0] seem most promising. 



Here we focus on the simplest solid-light model HI], de- 
scribing L optical microcavities each containing a single two- 
level atom with states ||), ||) separated by energy e. The 



JC 



P J 



(l) 



Hamiltonian reads = 1) 
H = - t J24a j + J2H:, 

Hf = elTiXTil+woaK+fldTiXiilOi + UiXTilol)- 

Here c^o is the cavity photon energy, and A = c^o — e defines 
the detuning. Each cavity is described by the well-known 
Jaynes-Cummings (JC) Hamiltonian H iC . The atom-photon 
coupling g (aj, a i are photon creation and annihilation op- 
erators) gives rise to formation of polaritons (combined atom- 
photon excitations) whose number N v = (a\a i J r\]i) 
is conserved and couples to the chemical potential \i [7] . We 
consider nearest-neighbor photon hopping with amplitude t, 
define the polariton density n = (N p )/L, use g as the unit of 
energy and set coo/g |23|], ks and the lattice constant to one. 

Hamiltonian (Q]) represents a generic model of strongly cor- 
related photons amenable to numerical methods. Existing 
theoretical work has focused on mean-field calculations fl|], 
exact diagonalization of few cavity systems y, [j, [ZD, an d 
the one-dimensional case |8|]. Here we employ a quantum 
many-body method for the thermodynamic limit to explore 
the physics of the model. In particular, quantum fluctuations 
on a finite length scale are included. We discuss Mott lobes 
(also at experimentally relevant finite temperatures), the effect 
of detuning and for the first time in such systems calculate 
single-particle spectra, a necessary connection to experiment 
and also a key metric in early proof of concept calculations of 
CAOL systems. 

The Variational Cluster Approach (VCA) — introduced first 
for strongly correlated electrons H 1 311 — has previously been 
applied to the BHM 11411 . The main idea is to approximate 
the self-energy 5] of the infinite system by that of a finite 
reference system. The matrix notation includes orbitals and 
bosonic Matsubara frequencies, T, = E a/ 3(icc; n ). The opti- 
mal choice for T, follows from a general variational principle 
Sft [S\ = 0, ft being the grand potential. Trial self-energies 
from isolated clusters (with L c sites) are parametrized by the 
one-particle parameters £ c = {t c , e c , cJq, fi c } of the reference- 
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system Hamiltonian, i.e. T, = £)(£ c ). For bosons |[l 



ft = ft c + Trln (G " 



Trln(G c ) 



(2) 



Here, ft c , G c , and X) are the grand potential, Green's func- 
tion, and self energy of an isolated cluster, and Go is the non- 
interacting (g = 0) Green's function. Go, G c , and T, are eval- 
uated at bosonic Matsubara frequencies iuj n , and traces in- 
clude frequency summation. The stationary solution is given 
by dfl/d£ c = 0. Traces can be evaluated exactly using only 
the poles of the Green's function but not their weig hts [13,141. 
The poles u m of Go -1 — 5] are obtained from a bosonic for- 
mulation of the Q-matrix method LI 5]. At temperature T > 0, 
the required matrix diagonalization restricts L c and T max . For 
simplicity, we restrict ourselves to a single variational param- 
eter £ C = 00 Q. 

The full quantum dynamics are taken into account exactly 
on the length scale of the cluster L c , and even for L c = 1 VCA 
results are beyond the mean-field solution . The present 

formulation cannot describe the properties of the SF phase, as 
the required symmetry-breaking term = if; (a\ + a i ) 
cannot be cast into a single-particle operator. However, this 
does not affect the accuracy of the phase boundary of the MI. 
The numerical effort is very moderate as compared to, e.g., 
the density matrix renormalization group [8], and the VCA 
provides T = and T > static and dynamic properties in 
one and two dimensions. 

The T = phase diagram in one (chain) and two dimen- 
sions (square lattice) is shown in Fig.Q] There exists a series 
of Mott lobes with integer polariton density n mt = 0, 1, . . . 
and compressibility k = dn/dfi = 0JLJ] - In contrast to re- 
cently proposed photonic MI phases 0], (a^a) fluctuates even 
for constant, integer n. Inside the lobes, where the VCA yields 
a solution, the system has an energy gap E g ^ (E g ,h) for adding 
a particle (hole) equal to the vertical distance of fi from the 
upper (lower) phase boundary 112411 . The spacing of the points 
ji n = ujo + g(V™ — V n + 1) where adjacent lobes touch at 
t/g = decreases quickly with increasing n, in contrast to 
the BHM where /i n = Un. Convergence with cluster size L c 
is surprisingly fast. The ID data agree well with exact results, 
although the VCA slightly underestimates t* (the value of t at 
the lobe tip, t*/g = 0.2 in (8D). FigureQJb) represents the first 
accurate (non-mean-field) phase diagram in 2D, arguably the 
most important case for experimental realizations. 

For t/g = 0, the MI states are EU^n);' where the n- 
polariton eigenstate of H ]C (the \ — ,n) branch in HI]) is a su- 
perposition of photonic (||, n), with n photons) and excitonic 
states (|t, n — 1)), 



|^n) = r](A, n) |I, n) + 0(A, n) ||, n-1) 



(3) 



Lobes with n > 1 are much smaller due to the effective polari- 
ton repulsion decreasing with n, and we focus on the n = 1 
case for which quantum effects are strongest, and which can 
be easily initialized experimentally. 

Figures [TJa) and (b) show significant deviations from the 
parabolic lobes predicted by mean-field theory in both one 
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FIG. 1: (Color online) Ground-state phase diagram: (a) ID, (b) 2D 
and (c) different coordination numbers z. Lines are guides to the eye. 



and two dimensions. However, quantum fluctuation effects 
diminish quickly with increasing coordination number z (Fig- 
ure QIc)). In particular, the reentrant behavior with increas- 
ing t/g [16] and the cusplike tip indicative of the Berezinskii- 
Kosterlitz-Thouless transition 111711 exist only for z = 2. 

Spectral properties play a key role in understanding con- 
densed matter systems, but are notoriously difficult to cal- 
culate accurately. In cavities, the occupation and spectra 
can be directly monitored through luminescence spectra, with 
angular emission translating into momentum k, much more 
straightforward than in CAOL. 

The VCA yields the single-particle spectral function 
A(k,u) = — 7r -1 ImG(/c, uj), shown in Fig. [2] for the pa- 
rameters marked in Fig. [T] Important for detection in ex- 
periment, the MI state is characterized by cosinelike parti- 
cle and hole bands, separated by the Mott gap E g , minimal 
at fe = 0, which decreases with increasing t/g, and even- 
tually closes in the SF phase. The phase boundaries fi± 
are related to the particle/hole dispersions £ p /h(& = 0, t) via 
= £ P (0,£) and fi-(t) = £h(0jX and the Mott gap is 
given by E g (t) = e p (0, t) - e h (0, t) El. 

For t <C t*, particles are much lighter than holes (e.g., 
tm p w 0.33, trrih ~ 0.81 inFig.[2tc); obtained from parabolic 
fits). The particle/hole bandwidth scales almost linearly with 



0.80 
0.60 
M 0.40 
i 0.20 



(d) 2D, I 



(e) 2D, II 



(f) 2D, III 



(0,0) (%,%) (ti,0) (0,0) (0,0) (%,%) (71,0) (0,0) (0,0) (7i,7i) (jt,0) (0,0) 
k k k 

FIG. 2: Single-particle spectra at T = in ID [(a)-(c), L c = 6] and 
2D [(d)-(f), L c = 2 x 2]. I-III refer to the marks in Fig.[l] i.e., (a) 
t/p = 0.16, (b) 0.08, (c) 0.01, (d) 0.05, (e) 0.02, (f) 0.005. 



t, confirming the weakly interacting Bose gas picture for 
the superfluid fraction of particles/holes doped into the MI 
|[l7ll . The hole bandwidth in both one and two dimensions 
is about zt. For the BHM model, the ratio of the bandwidths 
is Wp/Wh = 2 for the MI with n = 1 due to the effective par- 
ticle (hole) hopping t v = (n + l)t (t h = nt) HEl . For the 
JCM, the matrix elements for hopping of one particle or hole 
in the MI have to be evaluated using the dressed states |^ n ), 
and we find Wp/W h = 3/2 + y/2 « 2.91 in one dimension, 
in good agreement with Figs.[2a)-(c). The bandwidth ratio is 
fairly independent of t/U respectively t/g because the inter- 
action energy of a single particle or hole is the same at every 
site. In the BHM, emergent particle-hole symmetry leads to 
m$ ~ —> for t — » t* [20]. This behavior is also seen 
in the JCM (m p /m h « 1.3 for t = 0.16 in ID), but the VCA 
does not permit a detailed analysis near t* . 

In contrast to CAOL, the parameters e and c^o and hence 
the detuning A = e — c^o can easily be changed experimen- 
tally. This permits to tune the character of the bosonic parti- 
cles. The states \ip n ) are excitonlike for cj ^> photonlike 
for coo <C e, and polaritonlike for coo w e. The ID phase 
diagrams for A jg = ±2 are shown in Fig.[3a),(b). The exci- 
tonic system ( A/g = 2) with small photon-mediated hopping 
exhibits a large n = 1 Mott lobe, whereas the latter is very 
small in the photonic system with small exciton-mediated in- 
teraction. Reentrant behavior due to quantum fluctuations is 
seen for the excitonic case, but is absent in the photonic case 
for large enough L c . 

Single-particle spectra in one dimension near the lobe tips 
are shown in Fig.Oc),(d). The excitonic system shows a very 
large ratio of particle and hole bandwidths (Wp/Wh ~ 7), 
whereas Wp/Wh ~ 2.1 (very close to the BHM, since « 
||, 1)) for the photonic case. These values result from the dif- 
ferent admixture of the states ||, 0) and ||, 1) to |^i) depend- 
ing on A. In particular, the approximate relation Wh « zt 




FIG. 3: (Color online) Phase diagram for detuning (a) A/g — —2, 
(b) A/g = 2, and spectra for (c) A/g = 2, t/g = 0.7, and (d) 
A/g = -2, t/g = 0.0.018 in ID at T = 0. 



found for A = does not hold. The excitonic MI with 
n = 1 is approximately given by f[ i | j, 0) i9 whereas we have 
Y[ i ||, 1)^ for the photonic case. Finally, the incoherent fea- 
tures in Fig.Oc) may originate from the finite cluster size L c , 
and should be addressed using other exact methods. 

Solid state quantum devices will operate at finite tempera- 
ture, much higher than in CAOL (estimated at nK). This leads 
to the important question of the stability of the MI state at 
T > 0, which will ultimately determine their technological 
usefulness. Strictly speaking, there is no true MI at T > 
due to thermal fluctuations. However, there exist regions 
where fluctuations are small enough for the system to be- 
have like a MI for experimental purposes 112111 . We determine 
the region of existence of the MI using the stringent criterion 
An = \n — nint| < 10 -4 , corresponding to a "worst case sce- 
nario" since experimentally the MI will survive as long as the 
Mott gap ~ g is large compared to thermal fluctuations. 

Results for A = are shown for one and two dimensions 
in Fig. [4] Finite- size effects diminish quickly with increasing 
temperature (see Fig. |4^a)). The size of the n — 1 Mott lobe 
is significantly reduced with increasing temperature, and An 
exceeds 10 -4 at T* / g w 0.03 in both one and two dimen- 
sions. This value is consistent with the onset of deviations of 
n from 1 in the exact atomic-limit results, and agrees with the 
few-cavity results of (70. Lobes with n > 1 disappear at much 
lower temperatures (not shown). In contrast to T = 0, Mott 
lobes with n and n ± 1 do not touch at t/g = 0, but are sur- 
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FIG. 4: (Color online) Mott-like regions for different T. 

rounded by the normal fluid (NF) phase for sufficiently small 
t/g, and a transition NF-SF occurs at larger t/g. However, 
while the VCA can distinguish between the MI and the NF, 
the phase boundary to the SF at T > cannot be determined 
accurately. The gap in A(k : uj) increases with T 112 ill . Parti- 
cle/hole bands are still well defined at T > 0, and spectra (not 
shown) look very similar to T = 0. Excitation spectra of the 
NF and SF deserve detailed future studies. 

The quasi-MI lobes in Fig. [4] are dominated by the point 
(t = 0, ji « 0.3) where E g is maximal and by the T = 
lobe tip (/i* , t* ) . The MI is destroyed by thermal excitation of 
particles or holes, and therefore survives longest near fi = 0.3 
where the energy cost is largest. Besides, T* is determined 
by the size of the lobe at T = 0, t = 0. Hence, the system 
at T > is dominated by atomic-limit physics (yielding the 
same T* in ID and 2D) and the T = fixed point (quantum 
critical point) of the MI-SF transition [17]. For the ID BHM 
d, the VCA yields T*/U « 0.059, lower than T*/U = 0.2 
found in J22I1 using a different criterion ll25ll . The ratio of the 
critical temperatures for the ID JCM and ID BHM is 0.56, 
close to the ratio of the t = 0, T = Mott gaps (0.59). 

Solid state cavity-QED systems offer the possibility of large 
g. Taking a conservative estimate g = 10 10 Hz[l] gives for 
the effective temperature of the quantum system T* > 14 
mK. Even taken as the actual temperature, this is well ac- 
cessible experimentally. The hopping rate 1/t has to be fast 
enough to permit equilibration before photon loss and de- 
phasing set in. Indeed, taking t/g = 0.01 (inside the Mott 
lobe in Fig. HJb)), we obtain a realistic t~ x > 0.14 ps. Fi- 
nally, T* for the n = 1 MI is enhanced by detuning A > 
(Fig.Oa)), which additionally increases t* and |, estimated 
as (/A £*) « (0.08, 0.02) meV in 2D for T = A = 0. 



In summary, motivated by theoretical predictions and on- 
going experimental advances, we have studied polariton Mott 
phases by means of a versatile quantum many-body approach. 
Phase diagrams, single-particle spectra and finite tempera- 
ture effects have been related to the known features of Bose- 
Hubbard models and possible experiments, and the novel 
physics emerging from detuning has been explored. 

This work was supported by the FWF Schrodinger Fellow- 
ship No. J2583, the DFG research unit FOR538, and by NSF 
Award No. DMR-0502047. 



[1] A. D. Greentree, C. Tahan, J. H. Cole, and L. C. L. Hollenberg, 

Nat. Phys. 2, 856 (2006). 
[2] M. J. Hartmann, F. G. S. L. Brandao, and M. B. Plenio, Nat. 

Phys. 2, 849 (2006). 
[3] D. G Angelakis, M. F. Santos, and S. Bose, Phys. Rev. A 76, 

031805(R) (2007). 
[4] M. J. Hartmann and M. B. Plenio, Phys. Rev. Lett. 99, 103601 

(2007); Phys. Rev. Lett. 100, 070602 (2008). 
[5] N. Na, S. Utsunomiya, L. Tian, and Y. Yamamoto, Phys. Rev. 

A77,031803(R) (2008). 
[6] S.-C. Lei and R.-K. Lee, Phys. Rev. A 77, 033827 (2008). 
[7] M. I. Makin, J. H. Cole, C. Tahan, L. C. L. Hollenberg, and 

A. D. Greentree, Phys. Rev. A (accepted). 
[8] D. Rossini and R. Fazio, Phys. Rev. Lett. 99, 186401 (2007). 
[9] M. Lewenstein, A. Sanpera, V. Ahufinger, B. Damski, A. Sen, 

and U. Sen, Adv. Phys. 56, 243 (2007). 
[10] D. E. Chang, V. Gritsev, G. M origi, V. Vuletic, M. D. Lukin, 

and E. A. Dernier. larXiv:07 12. 1 8 1 71 (unpublished). 
[11] J. Koch, T. M. Yu, J. Gambetta, A. A. Houck, D. I. Schus- 
ter, J. Majer, A. Blais, M. H. Devoret, S. M. Girvin, and R. J. 

Schoelkopf, Phys. Rev. A 76, 042319 (2007). 
[12] C. W. Lai, N. Y. Kim, S. Utsunomiya, G Roumpos, H. Deng, 

M. D. Fraser, T. Byrnes, P. Recher, N. Kumada, T. Fujisawa, 

and Y. Yamamoto, Nature 450, 529 (2007). 
[13] M. Potthoff, M. Aichhorn, and C. Dahnken, Phys. Rev. Lett. 91, 

206402 (2003); M. Potthoff, Eur. Phys. J B 36, 335 (2003). 
[14] W. Roller and N. Dupuis, J. Phys.: Condens. Matter 18, 9525 

(2006). 

[15] M. Aichhorn, E. Arrigoni, M. Potthoff, and W. Hanke, Phys. 

Rev. B 74, 235117 (2006). 
[16] T. D. Kiihner and H. Monien, Phys. Rev. B 58, R14741 (1998). 
[17] M. P. A. Fisher, P. B. Weichman, G. Grinstein, and D. S. Fisher, 

Phys. Rev. B 40, 546 (1989). 
[18] N. Elstner and H. Monien, Phys. Rev. B 59, 12184 (1999). 
[19] D. van Oosten, P. van der Straten, and H. T. C. Stoof, Phys. Rev. 

A 63, 053601 (2001). 
[20] B. Capogrosso-Sansone, N. V. Prokof'ev, and B. V. Svistunov, 

Phys. Rev. B 75, 134302 (2007). 
[21] D. van Oosten, P. van der Straten, and H. T. C. Stoof, Phys. Rev. 

A 67, 033606 (2003). 
[22] F. Gerbier, Phys. Rev. Lett. 99, 120405 (2007). 
[23] As in |[l|]; 000 is not an absolute frequency, so that this choice 

does not affect results. 
[24] For t/g = and density n, a particle/hole corresponds to one 

site having occupation n + 1 (n — 1). 
[25] There are no well-defined plateaus in the density even below 

T*/U = 0.2 in HI. 



